Wigner function and Schrodinger equation in phase space representation 
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We discuss a family of quasi-distributions (s-ordered Wigner functions of Agarwal and Wolf) and 
its connection to the so called phase space representation of the Schrodinger equation. It turns out 
that although Wigner functions satisfy the Schrodinger equation in phase space they have completely 
different interpretation. 
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I. INTRODUCTION 

Since the pioneering work of Wigner 0, generalized 
phase-space techniques have found useful applications 
in various branches of physics H, ^ |^. The main 
idea of this approach is to represent the density opera- 
tor p as a function (quasi-distribution) over the classical 
phase space This function fully characterized the 

quantum state and enables one to express the quantum- 
mechanical expectations as averages of classical observ- 
ables over the classical phase space. Moreover, it is the 
Wigner function which is directly related to the measure- 
ment. Then quantum tomographic methods 0,13 enable 
one to reconstruct the quantum state from the experi- 
mental data. Recently, the Wigner function was also ap- 
plied to study quantum entanglement and related issues 
for continuous systems (see e.g. jESEo^)- 

The Wigner function is only one particular example 
of such a quasi-distribution. Especially important role 
(e.g. in quantum optics) is played by the family of func- 
tions introduced by Cahill and Glauber ^11.] containing as 
the Wigner function, the Glauber-Sudarshan P-function 
[T^ITsll . and the Husimi Q- function In this paper we 
analyze another lesser-known family introduced by Agar- 
wal and Wolf Actually, all these quasi-distributions 
correspond to the particular quantization procedure, that 
is, different ordering of q and p, or, equivalently, differ- 
ent ordering of creation and a annihilation operators, 
respectively. 

The procedure of representing quantum states by 
quasi-distributions in phase space is closely related to 
the phase space formulation of quantum mechanics based 
on the noncommutative product known as Moyal prod- 
uct Tg', or more generally as a star-product 18] (see 
1^19, 20] for the compact formulation of the standard quan- 
tum mechanics in terms of the Moyal product). 

There is another phase space representation of quan- 
tum mechanics based on the works of Torres- Vega and 
Frederic 0, |23] (see also 113] )■ In this approach the 
(pure) quantum state is represented by the wave func- 
tion '0(r)j where T represents a point in phase space. It 
turns out that V'(r) satisfies so called "Schrodinger equa- 
tion in phase space". The quantity lV'(r)l^ is, therefore, 
treated as a probability distribution in phase space. This 
procedure was applied to study simple quantum systems 
l^HEIH' In a recent paper Li, Wei and Lii ^] found a 



general method of solving "Schrodinger equation in phase 
space" . 

The aim of the present paper is to relate the stan- 
dard phase space approach based on quasi-distributions 
functions to that of Torres- Vega and Frederic (TF). We 
show that one can easily produce the whole family of 
"Schrodinger equation in phase space" which is closely 
related to the family of s-ordered Wigner functions Ws 
of Agarwal and Wolf ,15J, that is, Wg are particular so- 
lutions of this family of equations. Now, according to 
the standard approach Ws defines the quasi-distribution 
in phase space, whereas the TF approach implies that 
IV's]^, where ips = 2TThWs, is a (true) probability distri- 
bution, i.e. IV-s]^ > and J \il)s? dV = 1- It should be 
stressed, that "Schrodinger equation in phase space" has 
an infinite number of solutions. It is a price one pays 
for using V(r) instead of •^{q). Each particular solution 
gives rise to the particular phase space representation of 
ordinary wave function Lp{q) in position representation. 

The paper is organized as follows. In Sec. ^we present 
general approach to phase space representation of the 
wave function and following |23 we discuss the general 
solution for "Schrodinger equation in phase space" . Sec- 
tion mil introduces the whole s-family of equations to- 
gether with the family of solutions. Then, after recall- 
ing the formulae for star-products in Section Hvl we show 
that s-ordered Wigner functions do solve the family of 
Schrodinger equations. We end with some conclusions in 
Section El 



II. 



SCHRODINGER EQUATION IN PHASE 
SPACE 



There is no a unique way to represent a quantum state 
as a wave function -0 = ''/-'(r); where F represents a point 
in a classical phase {q,p). In the standard approach one 
usually uses a coordinate ^p{q) or momentum ^p(jp) repre- 
sentations, respectively. To pass from ip{q) to ^(F) one 
has to invent an integral transformation 



V^(F) = / K{T;q')^{q')dq' , 



(1) 



where K{T\ q') denotes the integral kernel. Functions 
V'(F) defined by the above formula form a proper sub- 
space TLk of the Hilbert space L^(R^). The unitarity of 



transformation L^(R) 9 Lp{q) — > V'(r) G Ha' requires 

K{T;q')K{r;q")dT = 5{q' -q") , (2) 



where denotes a measure on the phase space. Clearly, 
there is a huge freedom in choosing K. Performing the 
following "gauge transformation" 



(3) 



with /(r) being an arbitrary real function, one obtains a 
new kernel still satisfying ((SJ- 

In the literature there are several well known examples 
of such a transform. Perhaps the most famous is the 
Bargmann (or Bargmann-Segal) transform defined by |28l 

in 



KB{z;q) = vr ^^'^ exp 



V2. 



zq 



(4) 



where z is a complex number, i.e. one uses R'^ = C. The 
corresponding space TYb of entire functions V'(^) equipped 
with the following inner product 



(V'llV': 



'2;b 



ipl{z)ij)2{z)dii(z) , 



(5) 



where dii{z) = 7r~^e~^^^ d'^z, is known as the Bargmann- 
Segal representation of the Hilbert space. 

A closely related kernel is connected to the coherent 
states representation j30l | 



Kcs{q,P;q')=^{r\q') 



(6) 



where 



(ri-z') - 



exp 



2A2 



(7) 



with |r) denoting the standard Glauber coherent state 
corresponding to 



{\-^q + i\h-^p)/V2 



(8) 



The parameter A is a natural length scale defined by the 
mass and frequency of the oscillator, i.e. A = hj iiuj. 
The corresponding Hilbert space Hcs carries the follow- 
ing inner product 



with 



(V'i|V'2)cs = / v^i*(r)V2(r)dr , 

dqdp 



dV:^ 



2Trh 



(9) 



(10) 



From now on we shall use the dimensionless convention 
(|1(J|I for dT. Note, that in this convention tp{r) is dimen- 
sionless, whereas the kernel K{T;q') has the same dimen- 
sion as the wave function in the position representation 



Actually, the formula lO was a starting point in con- 
structing phase space representation ofguantum mechan- 
ics of Torres- Vega and Frederick |^ |22J . They showed 
that if (p{q) satisfies the standard Schrodinger equation 



^n-gl^iq^t) 



viq.t) , (11) 



then ip{T) obtained from (p{q) via Kcs{T;q) satisfies the 
following "Schrodinger equation in phase space" 



z?i^^(r,t) = 



2m oq^ 



zh— 
dp 



VXr,i) . 



(12) 



Actually, performing the gauge transformation 'ip'{q,p) 
g-*P9/2?i^^g p'j one finds more symmetric formula 



1 

2m 



^'(r,i) 



(13) 



where 



P = 



p ..d 
2 dq ' 



(14) 



satisfy [Q, P] — ih and, therefore, they define phase space 
representation of position and momentum. This particu- 
lar representation corresponds to the gauge transformed 
coherent states kernel e~'^P'^^^^Kcs{q,Pi q')- 

Recently, the following stationary Schrodinger equa- 
tion 



2m V2 dq 



9 I -fc ^ 



V [^ + ih^ 1 V'(r) = E^piT) , 

(15) 

was postulated in |23. Now, V'(r) denotes an arbitrary 
phase space representation, that is, the integral kernel 
K{T; q) in ^ is not specified. The general solution of 
l(T3|l reads as follows (eq. (11) in pTfl 



V-(r) = e-*«P/2'"' / g{y)ip{q + y)e-iPydy , 



(16) 



where g{y) is an arbitrary nonzero function and ip{q) is 
the eigenfunction of the Schrodinger equation in coordi- 
nate representation corresponding to the eigenvalue E. 
Note, that the function g{y) uniquely defines an integral 
kernel Kg by 

Kg{q,p; q') = e-^^^'^' -"/'^^'^ g{q' - q) . (17) 



Note, that 



K{T;q')K{T;q")dr 

^S{q'^q") f g*{q'~q)g{q"^q)dq, (18) 



2 



and hence unitarity condition ^ implies 

\g{y)\^dy^l. (19) 

In particular, the following Gaussian 

g(y) = (7rA2)-i/4,-.V2A^ ^ (20) 

does satisfy (|19|l and one finds for the corresponding 



III. A FAMILY OF THE SCHRODINGER 
EQUATION IN PHASE SPACE 

Let us observe that the representation p4|l may be 
generalized to the whole family of representations. It is 
convenient to scale phase space variables (<?,_p) — ^ {2q, 2p) 
and to introduce 



(21) 



with s e R. One can easily verify that [Qs,Ps\ = ih 
for s ±1. Note however, that lims^±i[Qs, Ps] = 
Therefore, the values s = ±1 will be understood as ap- 
propriate limits. 

In analogy to (|15|l let us postulate the following family 
of Schrodinger equations 



X Vs(r) = £;v^,(r) . 

To solve this equation assume that 

-2i pq 



Ah d 
^'-'hd-p 



ipsiq,p) = exp 



h 1 + s 

One obtains the following equation for 
sr _ / 2q 



(f'siq,?) ■ 



(22) 



(23) 



h (1 + S)2 



2m 



V 



1 



(1 



ih d 



X 0s(g,p) = E(j)s{q,p) 



(24) 



Now we expand the potential F as a Taylor's series about 
{ih{\ — s) /2)d / dp for given q and use the partial Fourier 
transform 



Uq,p)^ / x.(9,y)e-^^^/'' . 



(25) 



Further, multiplying both side of (|24l) by exp{— ipy'/S} 
and integrating over p, one obtains 



(1 + 5)^ 

2m 4 dq 

= Exs{q,y) , 



V 



2 1-s 

q-\ — 7;—y 



s + l 



Xsiq,y) 

(26) 



which defines the standard Schrodinger equation in the 
^-representation 



1 - s 



-q 



-y 



(27) 



s + l" 2 

Therefore, the general solution of (|26() reads as follows 

Xs{q.y) ^ 9s{y)^{£.) . (28) 

where gs — gs{y) is an arbitrary nonzero function and 
(p(^) satisfies 



v{0 



^(0 = E^{^) 



(29) 



2m de 

Finally, the general solution of H23|) has the following form 

2i pq 



ipsiq^p) = exp 



X 9s{y)v 



h s + l 



dye 



1-s 



-ipy/h 



-q 



-y 



,s + i^ ■ 2 

Clearly, for each gs{y) it defines a family of kernels Kg 
Mq,P) = J Kl{q,p-q')^{q')dq' 

given by 

Kl{q,p] q' 



(30) 



(31) 



1-s 
X exp 



9s 



2q' 



4q 



1-s 1 - s2 

2i p{q'^q) 



1 - s 



(32) 



Again, the requirement of unitarity ^ implies the fol- 
lowing condition for the function gs'. 



\gs{y)\^dy = 



1 + s 



(33) 



Note that for s = 1 the formula for the kernel consider- 
ably simplifies 

Kl=' (q, p;q') = e-'^'/'"' gip)Siq " q') , (34) 

where g stands for the Fourier transform of g. In this 
case the ave function ip{q) has the following phase space 
representation 



,Piq,p) = e~^P'^/'^gipMq) 



IV. STAR PRODUCT 



(35) 



Now we show that the family of equations H22|l is closely 
related to the family of quasidistributions function in 
phase space. 

To describe all quantum phenomena in phase space, 
we have to determine a relationship between operators 
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and functions on the classical phase space. This corre- 
spondence is of course not unique. The most famous is 
based on the Wigner-Weyl transform J-^/w- if f{Q:P) is a 
phase space function then one defines the corresponding 
operator F 



F — J^wwif) 



by 



F = da dTf{a,T)e 



i{a-q+Tp) 



(36) 



where / denotes the Fourier transform of / 

fir, a) = ^ J da J dr f{q,p)e-^'^^^+^^'^ . (37) 

The inverse transform, i.e. F ^ f /is defined as follows 



f{q,p) ^ I dy {q- -y 



F 



ipy/h 



(38) 



If F corresponds to a density operator p then its inverse 
Wigner-Weyl transform recovers celebrated Wigner func- 
tion W = T^\,{p) 



^il^P) = ~^ J dy {q - y\p\q + y)e 



2ipy/h 



(39) 



Now, the noncommutative multiplication of operators 
introduces the following noncommutative multiplication 
of functions 



Fi-F2 = Tww{h*f2) 



(40) 



The formula for the star-product was derived long ago 
by Groenewold and Moyal [TfiHT^: 

/i * /2 = /i exp I y {l)gdp - 'dp'dg^ I /2 , (41) 



where d { d ) act on the left (right) side. The Moyal 
product is associative 

(/i * /2) * /a = /i * (/2 * fs) , 
but it is noncommutative 

/i * /2 7^ /2 * /i • 

Recall, that the operator e""^^ acts as a generator of trans- 
lations: e""^^ f{x) = f{x + a). Therefore, the defining 
formula H41|l may be rewritten in the following form 



/i(«,p)*/2(g,p) = fi\q + -^dp,p- -T^dqj f2{q,p) ■ 

(42) 

Now we can make crucial observation: equation 122|) for 
s = has the following form 



(43) 



where H{q,p) is the classical Hamilton function. 

It turns out that the similar structure may be estab- 
lished also for s 0. Let us introduce the following family 
of Wigner-Weyl transforms: 



by 



F,^ I da I dTf{a,T)e'^''^+''P'^e-'"'^/^ . (44) 



Clearly, for s — one recovers (I36II . This formula enables 
one to introduce the family of star-products 

-4 -^ww(«*«^) > (45) 
which reduces to (|40|l for s = 0. One easily finds 

a-ks b ^ aexp ^{1 - s) d qd p - {1 + s) d pd b , 

(46) 

which is equivalent to 

a{q,p)-ks b{q,p) = 

q+{l-s)'^dp,p-{l + s)'-^d^b{q,p) . (47) 

Therefore, the family of Schrodinger equations in phase 
space may be rewritten as follows: 



H *s = Ei/^s 



(48) 



This shows that the family of equations (12211 which is a 
direct generalization of equations used in |2lj and is 
closely related to the noncommutative structure induced 
by the family of star-products. 



V. WIGNER FUNCTION VS. PHASE SPACE 
WAVE FUNCTION 

Both Wigner function W{q,p) and the wave function 
tpiT) are objects defined on the classical phase space. 
V'(r) satisfies (gg 



H-kil^ = Eil} 



(49) 



What about Wl It turns out that the stationary Wigner 
function is uniquely determined by the following two 
equations 19] 



HkW = Wi^H ^ EW 



(50) 



that is, W satisfies the same equations as ip and addi- 
tionally it fulfills W -k H = EW which is equivalent to 



1 

2m 



ih 



= EW{q,p) . 



Wiq,p) 



(51) 
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We stressed, that there are infinite solutions of H49|l and 
there is only one solution of (|50|l . Clearly, the Wigner 
function does belong to the solutions of Indeed, 
taking 



g(y) = ^*(-y/2) 
the formula (|16|l implies 



V(r) 



(9-f)^(^ 



(52) 



(53) 



that is, V'(r) = 2TThW, where W{q,p) is a Wigner func- 
tion corresponding to ip. Note, that g{y) defined via H52|l 
satisfies and hence the corresponding kernel does in- 
deed satisfy 

As an example let us compare the solutions of (|49|l 
and (|50|l for the harmonic oscillator. Taking g{y) as in 
ipn|l one obtains from lO [H IH IH the following 
formulae corresponding to nth energy eigenstate ipnil)'- 



with r given by ||SJ), whereas 

W^{q,p) = ^1^L„ (^) exp f-|^) , (55) 

where L„ denotes nth Laguerre polynomial. Due to 
|rp — H/huj, one has for the probability distribution 
of transition from (pn to the coherent state |r) 



l^n(r)p = l )"expf-|- 



(56) 



Clearly, both iV'nP and Wn depends only upon the os- 
cillator energy H and both are normalized according to 



(57) 



Moreover, it its easy to show 



(58) 



Now, the family of Wigner- Weyl transforms -^^vvw 
ables one to introduce the following family of Wigner 
functions: Wg = (•^ww)~^(p)' where p stands for the 
density operator. One finds 



X {q-{l 



ipy/h 



(1 



(59) 



which reduces to W{q,p) for s — 0. The family Ws was 
introduced by Cahill, Glauber, Agarwal and Wolf [lllllSt . 
It satisfies two basic properties: it is normalized 



J Ws{q,p)dqdp = 1 , 



and for any quantum observable F 

Tr{Fp) = / Ws{q,p)f^s{q,p)dqdp 



where /_s = (-^ww) ^i-^)- ^'^^ s = the last formula 
reproduces well known property of the Wigner function 

Tr(Fp) = J Wiq,p)f{q,p)dqdp . 

Moreover, Wg provides correct quantum marginals: 

dq Ws{q,p) = {p\p\p) 

dpWs{q,p) = {q\p\q) . (60) 



It turns out that stationary s- Wigner functions Ws are 
uniquely determined by 



(61) 



Equation (|48|l for ips has infinite number of solutions 
whereas the set of two equations (|61|l has only one so- 
lution 

V'.(r) = Jdy e-'py/'' 

X f* (^Q - '^-^y^ ^ (q + ^-^yj , (62) 

i.e. V's(r) = 27rfi.Ws((?,p). Therefore, Wg is only one 
particular solution of H48|) . It is easy to see that taking 



9s{y) = 'P* ( -^4—^ 



(63) 



in H3U|I one obtains ips(r) given by (|62|) . In particular for 
s = 1 one obtains so called Kirkwood-Rihaczek function 
^il^p) = W^s=i(9:?') which in the case of pure state (p 
reduces to 



Kiq,p)^e^P^^'^p*ip)piq) 



(64) 



It was introduced by Kirkwood [2l| as an alternative for 
the Wigner function. Then, in 1968, the same formula 
was rediscovered by Rihaczek [s^ in the context of sig- 
nal time-frequency distributions (see 33] for a useful re- 
view). Recently, this function was analyzed and applied 
in various contexts in [sl I3I IsIIstL IM Isgf . 



VI. DISCUSSION 

Both the phase space wave function ip(T) and s-ordered 
Wigner function Wg encode the entire information about 
the quantum state 'p{q). Due to the basic property 



(65) 



5 



some authors call |V'(r)p a probability distribution in 
phase space. Clearly, quantum mechanics does not al- 
low for a genuine probability distribution in q and p\ To 
interpret jV'lr)!^ correctly note that formula ^ may be 
rewritten as the following inner product 

^(r) = {^r\^) , (66) 

where 

M<l')^K*{r;q') . (67) 

Let us consider kernels defined by (|17|l . Then, due to Ijiyil . 
ifriQ') is a normalized wave function in the position repre- 
sentation. Therefore, ['(/'(r)^ is the probability density of 
transition from the state ip to state fr- In particular for 
the coherent state kernel Kcs one has |V'(r)p — |(r|(y9)p 
which defines the Husimi function for the state ip. 

Now, the Wigner function defined quasi-distribution 
such that / W dqdp = 1 . Since is a special solution of 
the Schrodinger equation in phase space one has 

W^{q,p) = -^{^r\v^) , (68) 
Zirn 

where is the Wigner function corresponding to ip. 



Formulae and imply 

ipriq') - 2ip*{2q - q') 6-2*^(9'-'?)/'^ . (69) 

It should be stressed that the phase formulation based on 
the wave function is restricted to pure states only 

whereas the approach based on Wigner function works 
perfectly for general mixed states p. Therefore, this ap- 
proach is much more general. Note, that for mixed states 
one has 

and the equality holds for pure states only. Therefore, for 
general mixed states {2Trh)'^ cannot be interpreted as 
a probability distribution. 
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